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Abstract 

We consider the problem of reconstruction of a non-linear finite-parametric model 
M = M p {x), with p = (pi, . . . ,p r ) a set of parameters, from a set of measurements 
fij(M). In this paper fj,j(M) are always the moments rrij(M) = j x 3 M p (x)dx. 
This problem is a central one in Signal Processing, Statistics, and in many other 
applications. 

We concentrate on a direct (and somewhat "naive") approach to the above 
problem: we simply substitute the model function M p (x) into the measurements 
Hj and compute explicitly the resulting "symbolic" expressions of fij (M p ) in terms 
of the parameters p. Equating these "symbolic" expressions to the actual mea- 
surement results, we produce a system of nonlinear equations on the parameters 
p, which we consequently try to solve. 

The aim of this paper is to review some recent results (mostly of [T3| \TS[ 
[T9| l28| l29l [30l I3T1 [47] ) in this direction, stressing the algebraic structure of the 
arising systems and mathematical tools required for their solutions. 

In particular, we discuss the relation of the reconstruction problem above with 
the recent results of [H El [Ml E3 E21 E3 S3] on the vanishing problem of general- 
ized polynomial moments and on the Cauchy-type integrals of algebraic functions. 

The accompanying paper [21] (this volume) provides a solution method for 
a wide class of reconstruction problems as above, based on the study of linear 
differential equations with rational coefficient, which are satisfied by the moment 
generating function of the problem. 
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1 Introduction 



In this paper we consider the following problem: let a finite-parametric family 
of functions M = M p (x), x G M. m be given, with p = (pi, . . . ,p r ) a set of 
parameters. We call M p (x) a model, and usually we assume that it depends 
on some of its parameters in a non-linear way (this is always the case with 
the "geometric" parameters representing the shape and the position of the 
model). 

The problem is: 

How to reconstruct in a robust and efficient way the parameters p from a set 
of "measurements" /ii(M), . . . , fii(M) ? 

In this paper \ij will be the moments rrij(M) = J x^ M p (x)dx. This as- 
sumption is not too restrictive - see, for example, [UH ITT] . 

The above problem is certainly among the central ones in Signal Process- 
ing (Non-linear matching), Statistics (Non-linear Regression), and in many 
other applications. See [III EH HH1 ESI ESI EDI ED HZ] and references there. 

We concentrate in the present paper on a direct (and somewhat "naive" ) 
approach to the above problem: we simply substitute the model function 
M p (x) into the measurements and compute explicitly the resulting "sym- 
bolic" expressions of Hj(M p ) in terms of the parameters p. Equating these 
"symbolic" expressions to the actual measurement results, we produce a sys- 
tem of nonlinear equations on the parameters p which we try to solve. 

Certainly, the polynomial moments do not present the best choice of 
measurements for practical applications since the monomials x^ are far away 
from being orthogonal (see, for example, [46J). However, the main features 
of the arising non-linear systems remain the same for a much wider class of 
measurements, while their structure is much more transparent for moments. 

The aim of this paper is to review some recent results (mostly of [TTJ E2 
IT8| IT9l |28| |29| 1301 |3T| ITT] ) in this direction, stressing the algebraic structure 
of the arising systems and mathematical tools required for their solutions. 
In particular, we stress the role of the moment generating function. 

We start with some initial examples of the models M p (x) in one dimen- 
sion: these are polynomials and rational functions. Then we consider linear 
combinations of (^-functions. The system which appears in this example is 
typical in many application. We discuss one of the classical solutions meth- 
ods, following [371 m\ E3 HE1 1211 w\. 

Next we deal with piecewise-solution of linear differential equations, pro- 
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viding some pre-requisites for the reconstruction method described in [21] 
(this volume). Then we consider piecewise-algebraic functions of one vari- 
able. We prove injectivity of the finite moment transform on such functions, 
and discuss the relation of the reconstruction problem for such functions 
with the recent results of [U El El EU ESI E2] on the vanishing problem of 
generalized polynomial moments. 

In two dimensions we shortly present results of [321 [TS], ED] concerning 
reconstruction of polygons from their complex moments, as well as results 
of [19] on reconstruction of "quadrature domains" . Finally we consider the 
problem of reconstruction of ^-functions along algebraic curves, relating it to 
the vanishing problem of double moments ([U EJ [221 [221 El 155]). 

We almost do not touch the classical Moment Theory, refereeing the 
reader to [33] and especially to [TSl [2QI ES HQl EH HS], where, in par- 
ticular, a review of the classical results and methods is given, as applied to 
the effective reconstruction problem. 

We also don't discuss in this paper the problem of noise resistance. It is 
treated in [EH CEH1 [2E1 ESI EI] ■ 

1.1 Applicability of the "direct substitution" method 

The key condition for applicability of our approach is the assumption that 
the signals we work with can be faithfully approximates by a priori known 
"simple" geometric models. 

A natural question is: to what extent this assumption is realistic? The 
answer to this question is twofold: 

1. In many specific application the form of the signal is indeed known a 
priori. Besides the wide circle of applications mentioned in [TH [T21 [TBI [2HJ 
|2"9"| EDI ED HZ] notice that this is usually the case in visual quality inspection. 
Similar situations arise in some medical applications where a non-linear para- 
metric model of an important pattern has to be matched to the radiology or 
ultrasound measurements data. 

2. A general applicability of our approach in problems involving image 
acquisition, analysis and processing depends on a possibility to represent 
general images of the real world via geometric models. 

The importance of such a representation in many imaging problems, from 
still and video-compression to visual search and pattern detection is well- 
recognized. Some initial implementations of geometric image "modelization" 
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have been suggested, in particular, in [27J El El E5] • See [T5] and references 
there for a general overview and analysis of the performance of edges-based 
methods in images representation. 

However, in general the "geometric" methods, as for today, suffer from 
an inability to achieve a full visual quality for high resolution photo-realistic 
images of the real world. In fact, the mere possibility of a faithful capturing 
such images with geometric models presents one of important open problems 
in Image Processing, sometimes called "the vectorization problem". 

Certainly, this current status of affairs makes problematic immediate 
practical applications of general imaging methods based on geometric model. 

Let us express our strong belief that a full visual quality geometric- 
model representation of high resolution photo-realistic images is possible. 
As achieved, it promises a major advance in image compression, capturing, 
and processing, in particular, via the approach of the present paper. 

Recently some "semi-linear" approaches have emerged providing a reli- 
able reconstruction of "simple" (and not necessarily regular) signals from a 
small number of measurements. In these approaches (see [HI [TO] and refer- 
ences there) "simplicity" or "compressibility" of a function is understood as a 
possibility for its accurate sparse representation in a certain (wavelet) basis. 

A somewhat more general approach to the notion of a complexity of func- 
tions has been suggested in [HHH]: here we take as a complexity measure the 
rate of semi-algebraic approximation. If the wavelet base is semi-algebraic, 
"compressible" functions have low semi-algebraic complexity. The same is 
typically true for functions allowing for a fast approximation by various types 
of non-linear models. 

2 Examples of moment inversion: one vari- 
able 

In this section we consider some natural examples of the models M p (x) in one 
dimension and of their reconstruction from the moments. These are poly- 
nomials, rational functions, linear combinations of ^-functions, and the class 
Ad of piecewise- analytic functions, each piece satisfying a fixed linear differ- 
ential operator D with rational coefficients. (Piecewise-polynomials belong 
to Ad for D = -fir). Then we consider piecewise- algebraic functions. 

In this paper we use as one of the main tools in solving the moment 
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inversion problem the moment generating function I g {z) defined as 



2.1 Vetterli's approach 

In jl?! EHl E] an important class of signals has been introduced, possessing 
a "finite rate of innovation", i.e. a finite number of degrees of freedom per 
unit of time. Usually such signals are not band-limited, so classical sampling 
theory does not enable a perfect reconstruction of signals of this type. In 
[4"Tl [28j [11] it was shown that using an adequate sampling kernel and a 
sampling rate greater or equal to the rate of innovation, it is possible to 
reconstruct such signals uniquely. The behavior of the reconstruction in the 
presence of noise has been also investigated. 

The main type of signals for which explicit reconstruction schemes have 
been proposed include linear combinations of ^-functions and their deriva- 
tives, splines, and piecewise polynomials. In spite of a somewhat different 
setting of the problem, the reconstruction schemes turn out to be mathe- 
matically similar to the ones presented below. In fact, moments enter, as 
an intermediate step, the reconstruction procedure in [11], and systems very 
similar to (2.7) and (2.9) below explicitly appear in [27J [281 HI]- ft is a re- 
markable fact (although traced at least to [37J) that exactly the same systems 
arise in exponential approximation ([IT]), in reconstruction of plane polygons 
( [131 ESI ED] , see Section 3.1 below), in reconstruction of quadrature domains 
([19J, see Section 3.2 below), in Pade approximations, and in many other 
problems. 

In [29] the approach of [UJ, [281 EH] is extended to some classes of paramet- 
ric non-bandlimited two-dimensional signals. This includes linear combina- 
tions of 2D ^-functions, lines, and polygons. Notice that the first problem, in 
its complex setting (where we consider as the allowed measurements only the 
complex moments //&(/) = f f z k f(x,y)dxdy) leads once more to a complex 
system (2.7). 




(2.1) 
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2.2 Polynomials 

Let P(x) be a polynomial of degree d, P(x) = Y^j=o a j x ^ ■ F° r the k-th 
moment rrik{P) we have 

n\ d d d 

m k {P) = / a ^ +kdx = E TTTTT = E ^ ' ( 2 ' 2 ) 

if we put hkj = ■ Now let a denote the column- vector of the coefficients 
dj of the polynomial P(x) and let m denote the column- vector of the moments 
m (P), . . . , rrid(P). We get the following linear system: 

Ha = m, H = (h kj ). (2.3) 

Notice that the matrix H is a Hankel matrix: the rows of this matrix are 
obtained by the shifts of its first row. More specifically, the matrix H belongs 
to the class of Hilbert-type matrices (see [ZE]). In particular, its determinant 
is nonzero, and system (2.2) has unique solution. Therefore, we have 

Proposition 2.1 A polynomial P(x) of degree d can be uniquely recon- 
structed from its first d+ 1 moments m (P), . . . , md(P), via solving system 
(2.3). 

Notice, however, that the smallest eigenvalue X m i n (H) behaves asymptoti- 
cally for d — > oo as follows: 

X mm (H) = KVdp-^ d+1 \l + o(l)), 

where K = 87rv / 27r24 and p = 1 + y/2. ([2T]). Therefore, the inversion of the 
matrix H becomes problematic for large d. 

Notice also that for each fixed polynomial P(x) expression (2.1) defines 
rrik(P) as a rational function of k. 

As for the moment generating functions, we have 

Proposition 2.2 Ip{z) = — - log(l — -) + P(-), with P(s) a polynomial of 
degree d — 1 in s. 

Proof: We have P(t) = P(t)(t - l ) + Pi 1 ) where P(t) is a polynomial of 
degree d — 1 in t with the coefficients - polynomials of degree d — 1 in -. 
Hence 

Integrating from to 1 now provides the required expression. 
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2.3 Rational functions 

Let R(x) be a rational function of degree d, R(x) = ^M, deg Q = d, deg P < 
d — 1 (we assume that R does not have a "polynomial part"). Thus 

d-i d 

p(x) = yja/a^, <5(x) = fejx- 7 . 

3=0 j=0 



We have 



Hence 



d 

P(x) = Q(x)R(x) = ^b j x j R(x) 



d 

m k{ P ) = ^2 b 3 m k+j{R), k = 0, 1, ... , 

3=0 

and using our notations from Section 2.1 above we finally get a system for 
the unknowns aj, bj 

d-l d 

hk i a i = Yl m k+j(R)bj, k = 0, 1, . . . , 2d, (2.4) 

j=0 j=0 

where, as above, hkj = We do not analyze here the solvability con- 

ditions for (2.4) (compare, however, Lemma 2 in [21]). Let us notice also 
that the counting of the sign changes as in Section 2.6 below shows that a 
rational function R(x) of degree d can be uniquely reconstructed from its 
first 4d moments m (R), . . . , m id (R). 

To compute the moment generating function Ir(z) let us assume that the 
roots oti, . . . , «rf of Q are all distinct. Then R(t) = Ylt=i ~i~a~ anc ^ denoting 

- by w we get t _ w = J2i=i (t- ai )lt-w) = Si=i ^(( ai - w )(t- ai ) ~ ( ai -w)(t-w)^ ■ 
Transforming integral (2.1) as in the proof of Proposition 2.2, and integrating 
we finally get 

Proposition 2.3 The moment generating function Ir(z) of a rational func- 
tion R(x) is given by 



i fl co = ^x;^^M— )], «> = -■ (2-5) 

/ — ' a,- —wou w z 



i=l 
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2.4 Linear combination of ^-functions 

Let g(x) = Tif =1 Ai5(x — Xi). For this function we have 

m k (g)= [ x k Y% =1 A i 5(x-x i )dx = Y% =l A i x k i . (2.6) 
Jo 

So assuming that we know the moments rrik{g) — k — 0, 1, . . . , 2n — 1, we 
obtain the following system of equations for the parameters A, and x^, i = 
1, . . . , n, of the function g: 

£™ =1 A4 = ct k , k = 0, 1, . . . , 2n - 1. (2.7) 

Notice that system (2.7) is linear with respect to the parameters Ai and 
non-linear with respect to the parameters Xi. 

System (2.7) appears in many mathematical and applied problems. First 
of all, if we want to approximate a given function f(x) by an exponential 
sum 

f(x) ss de aix + C 2 e a2X + ■■■ + C n e a " x , 

then the coefficients Cj and the values /u, = e ai satisfy a system of the form 
(2.7) with the right-hand side (the "measurements") being the values of f(x) 
at the integer points x = 1,2,.... (see [17], Section 4.9). The method of 
solution of (2.7) which we give below, is usually called Prony's method ([37]). 

On the other hand, system (2.7), recurrence (2.9) and system (2.10) below 
form one of the central objects in Pade approximation: see, in particular, [53] 
and references there. 

System (2.7) appears also in error correction codes, in array processing 
(estimating the direction of signal arrival) and in other applications in Signal 
Processing (see, for example, [30], [TTJ and references there). 

In [131 [T5] [30] system (2.7) appears in reconstruction of plane polygons 
from their complex moments. These results are shortly described in Section 
3.2 below. 

This system appears also in some perturbation problems in nonlinear 
model estimation. 

We give now a sketch of the proof of solvability of (2.7) and of the solution 
method, which is, essentially, the Prony's one. We follow the lines of [30, 25J. 
See also a literature on Pade approximation, in particular, [33] and references 
there. 
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Theorem 2.1 A linear combination g(x) of n 5-functions can be uniquely 
reconstructed from its first 2n — 1 moments mo(g), • • • , m2 n -x{g), via solving 
system (2.7). 

Proof: Representation (2.6) of the moments immediately implies the follow- 
ing result for the moments generating function I g (z): 

Proposition 2.4 For g(x) = Y^ =1 Aib~(x — Xj), the moments generating func- 
tion I g (z) is a rational function with the poles at Xi and with the residues at 
these poles A^: 

'« = E "-r^r < 2 - 8 > 

We see that the function I(z) encodes the solution of system (2.7). So to solve 
this system it remains to find explicitly the rational function I(z) from the 
first 2n its Taylor coefficients oeo, . . . , ctm-x- This is, essentially, the problem 
of Pade approximation (|33j). 

Now we use the fact that the Taylor coefficients of a rational function of 
degree n satisfy a linear recurrence relation of the form 

m r+n = E^Cjm r+ j, r = 0, 1, . . . . (2.9) 

Since we know the first In Taylor coefficients a^, . . . , a2„_i, we can write a 
lin ear system on the unknown recursion coefficients Cf. 

T^Cw+r = a n+r , r = 0, 1, . . . , n - 1. (2.10) 

Solving linear system (2.10) with respect to the recurrence coefficients Cj 
we find them explicitly. For a solvability of (2.10) see [171 [30l [25]. Now 
the recurrence relation (2.10) with known coefficients Ci and known initial 
moments allows us to easily reconstruct the generating function I g (z) and 
hence to solve (2.7). 

Remark. Another proof of Theorem 2.1 can be obtained in lines of the proof 
of Theorem 2.2 below. Indeed, a difference of two linear combinations of n 
(^-functions can have at most 2n — 1 "sign changes" . Then we apply Lemma 
2.2. 
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2.5 Piecewise-solutions of linear ODE's 

In this paper we do not consider separately the case of piecewise-polynomials. 
See [37J where a method for reconstruction of piecewise-polynomials from 
samplings is suggested (which starts with a reconstruction of linear combi- 
nations of ^-functions and of their derivatives). Instead we consider, as a 
natural generalization of piecewise-polynomials, the class Ad of piecewise- 
analytic functions, each piece satisfying a fixed linear differential operator 
D with rational coefficients. Such functions are usually called "L-splines" 
(see |44[ |4"5] and references there). For piecewise-polynomials of degree d we 
have D = ■ Notice that Vetterli's method ([31]) can be extended also to 
our class Ad- However, in the present paper we stress another approach to 
the moment reconstruction problem for the class Ad- It is presented in the 
accompanying paper [23] (this volume), while here we provide a necessary 
background. 

Consider the equation 

Dy = y {k) + a k ^{x)y {k - l) + ■■■ + ai {x)y' + a {x)y = (2.11) 

with the coefficients dk-i(x), . . . , ao(x) real-analytic and regular on [0,1]. 
All the solutions of (2.11) on [0, 1] form a linear space Ld with the basis 
yi(x), . . . , yk{x) being the fundamental set of solutions of (2.11). For D = 
the space Ld consists of all the polynomials of degree at most k — 1, and we 
can take {yi(x), . . -,y k (x)} = {l,x,x 2 , . . . ,x k ^}. 

Now we consider the class Ad of all the piecewise-continuous functions 
g(x) on [0, 1] with the jumps at x\, . . . ,x n G [0, 1], such that on each conti- 
nuity interval Aj = the function g(x) satisfies Dg = 0. We extend 
g(x) by the identical zero outside the interval [0, 1]. 

We can represent g(x) on the intervals Aj in a "polynomial form": g(x) = 
Y^!j=i a ijyj{ x )i where yi(x), . . . ,yk(x) is the fundamental set of solutions of 
(2.11). Alternatively, we can parametrize g(x) on the intervals Aj by its 
initial data at the point Xj. We can further define "splines" of a prescribed 
smoothness in Ad- The constructions of [37J can be extended to this case. 

While till this point we could restrict our presentation to the real domain, 
in what follows it will be necessary to extend the consideration to the complex 
plane. 
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First we recall shortly some classical facts related to the structure of linear 
differential equations in the complex domain (see, for example, [221 H3] where 
these fact are presented in a form convenient for our applications). 

Consider the equation 

Dy = + a k - 1 (x)y«'- 1 ') + ■■■ + ai {x)y' + a {x)y = (2.12) 

with the coefficients a,k-i(x), . . . , ao(x) regular and univalued in the complex 
domain Q = C \ {x , . . . , x m }. We do not specify at this stage the character 
of possible singularities of at the points x , . . . , x m . 

The following proposition (see, for example, [13]) characterizes multival- 
ued analytic functions which are solutions of a certain equation of the form 
(2.12): 

Proposition 2.5 Any solution y(x) of (2.12) is a regular multivalued func- 
tion in Q, satisfying the following additional property (F): For any point 
w G Q the linear subspace L w spanned by all the branches of y(x) at w in the 
space 0(w) of all the analytic germs at w, has dimension at most k. 

Any regular multivalued function v(x) in Q with the property (F) satis- 
fies a certain equation of the form (2.12) of order at most k with all the 
coefficients regular and univalued in the domain Q. 

Let us remind that for a given function g(x) on [0, 1] the moment generating 
function I g (z) = Y^T=o m k{9) zk * s given by the Cauchy-type integral 

f 1 gifydt fgMdt 1 . 

W)=/ i 7 = w 7' w = ~- 2 - 13 

J 1 - zt J w - t z 

Now one of the basic classical facts about Cauchy-type integrals is that if g 
(on each its continuity interval) satisfies a certain equation of the form (2.12) 
then I g (z) satisfies another equation of this form. A proof (in a specific case 
which we need in the present paper) can be found in [361 H3] • bn these papers 
also specific ramification properties of I g (z) are studied for g algebraic. 

Now, in the accompanying paper [23] (this volume) the functions g(x) 
from the class Ad are considered. A non-homogeneous equation of the form 
(2.12) for I g (z) is presented explicitly, and on this base a reconstruction 
procedure is suggested. 
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2.6 Piecewise-algebraic functions 

Exact reconstruction of piecewise-algebraic (= semi-algebraic) functions can 
be considered as one of the ultimate goals of our approach. If we extend 
this class SA to SA(i/ji, . . . ,ipi), adding a finite number of fixed "models" 
x/ji, . . . ,tpi and allowing for all the elementary operations and for solving 
equations, we shall probably cover all the examples of interest. In partic- 
ular, such extensions include linear combinations of shifts and dilations of 
ipi , . . . , ipi - an important class appearing in reconstruction of signals with 
finite innovation rate f [47[ 1281 [TT]). and in wavelets theory. Extensions of 
this sort are also closely related to what appears in theory of o-minimal 
structures - see, for example, [12]. Because of the "finiteness results" in this 
theory we can hope that the "finite moments determinacy" of semi-algebraic 
functions (Theorem 2.2 below) can be extended to at least some important 
classes SA(ipi, . . . 

Let us remind that g(x) is an algebraic function (as usual, restricted to 
[0, 1]) if y = g(x) satisfies an equation 

a n (x)y n + a n _i(x)y n_1 + • • • + a x {x)y + a (x) = 0, (2.14) 

where a n (x), . . . , a (x) are polynomials in x of degree m. d = m + n is, by 
definition, the degree degg of g. 

We shall need the following simple properties of algebraic functions: 

1. The number of zeroes of an algebraic function g(x) defined by (2.14) 
does not exceed m (and so it does not exceed its degree deg g = m + n). 

2. A sum g(x) = gi(x) + g2^x) of two algebraic functions of degrees d\ 
and c?2 is an algebraic function, with the degree degg < 77(^1, g^)- 

We consider piecewise-algebraic functions on [0, 1]. Let such a function 
g(x) be represented by the algebraic functions g q {x) of the degrees d q , re- 
spectively, on the intervals A q = [x q ,x q+ i], q = 0,...,r, of the partition of 
[0, 1] by xo = < x\ < ■ • • < x r < x r+ i = 1. We define the combinatorial 
complexity, (or the degree) o~(g) of g as follows: 

Definition 2.1 (See J^g| \4^j)- The combinatorial complexity a(g) is the 
sum Yfq=i d Q + r - 

The specific choice of this expression is motivated by the following simple 
observation: the number of sign changes of a piecewise-algebraic function g 
on [0, 1] does not exceed o~(g). This follows directly from property (1) above. 
We need also the following lemma: 
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Lemma 2.1 Let g±,g 2 be piecewise- algebraic functions with a(gj) < d, j = 
1,2. Then for g = g\ ± g 2 the combinatorial complexity a(g) satisfies a{g) < 
n(d) = 2d(i](d, d) + 1), where i](d, d) is given by property (2) above. 

Proof: g has at most 2d jumps, and on each continuity interval its degree 
is bounded by rj(d, d). 

Now we can show that piecewise-algebraic functions are uniquely defined 
by their few moments. We do not touch in this stage the question of how such 
a function can be actually reconstructed from the moments data, postponing 
this problem till Section 2.6.1. 

Theorem 2.2 A piecewise-algebraic function of a combinatorial complexity 
d is uniquely defined by its first n(d) moments. 

Proof: Assume, in contrary to the statement of the theorem, that there are 
functions g\ and g 2 of complexity at most d, with exactly the same moments 
up to order s = n(d). Hence for the difference g = g 2 — gi ^ we have the 
vanishing of the moments up to s: m,j(g) = 0, j — 0, 1, . . . , s. By Lemma 
2.1 we have for the combinatorial complexity of g the bound a(g) < s. 
Consequently, the number of sign changes of g does not exceed s. The next 
trick comes from the classical moment theory: 

Lemma 2.2 // the number of the sign changes and zeroes of g(x) ^ does 
not exceed s then some of its first s moments m,j(g), j = 0, 1, . . . , s do not 
vanish. 

Proof: We can assume that g changes its sign at certain points t±, . . . ,ti, I < 
s, and preserves the sine between these points. Let us construct a poly- 
nomial Q(t) of degree / with exactly the same sign pattern as g: Q(t) — 
±(x — ti)(x — t 2 )---(x — ti). Write Q as Q(x) = Yli a j xj - We have 
g(x)Q(x) > everywhere, besides ti,...,ti and possibly some other iso- 
lated points. Therefore J* 1 g(x)Q(x) > 0. On the other hand, this integral 
can be expressed as a linear combination of the moments: g(x)Q(x) = 

ctj Jq xig(x)dx = Yli a j m j(9)- Hence some of the moments of g up to 
/ < s-th do not vanish. This proves Lemma 2.2. To complete the proof of 
Theorem 2.2 it remains to notice that the difference g = g 2 — g\ ^ on at 
least one of its continuity intervals. 
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2.6.1 Explicit moment inversion for algebraic functions 



As far as an explicit inversion of the moment transform of algebraic functions 
is concerned, we are not aware of any general approach to this problem. 
Piecewise-algebraic functions belong to the class Ad, as defined in Section 2.3 
above. However, the problem is that we do not know a priori the differential 
operator D which annihilates a given algebraic function g. (The form of D is 
known, but not the coefficients of the rational entries of D). This fact seems 
to prevent a direct application of the method of [24j to piecewise-algebraic 
functions. 

Let us analyze in more detail one special case. Assume that the algebraic 
curve y = g(x) is a rational one. This means that it allows for a rational 
parametrization 



The moments Trik(g) given by rrik(g) = L x k g{x)dx, k = 0, 1, ... , now can 
be expressed as 



where p denotes the derivative of P and = P(a), 1 = P(b). Moments 
of this form naturally appear in a relation with some classical problems in 
Qualitative Theory of ODE's - see [21 H El E] , [33]- EE]. 

Our problem can be reformulated now as the problem of explicitly finding 
P and Q from knowing a certain number of the moments in (2.16). 

Of course, in general we cannot expect this system of nonlinear equations 
to have a unique solution. Indeed, while the function y = g(x) is determined 
by its moments in a unique way, the rational parametrization P, Q of this 
curve in general is not unique. In particular, let W(t) be a rational function 
satisfying W(0) = 0, W{1) = 1. Substituting W{t) into P and Q we get 
another rational parametrization of our curve: 



x = P(t), y = Q(t), with P(t) = P(W(t)), Q(t) = Q(W(t)). (2.17) 



Consequently, we can ask the following question: 

Are all the solutions of (2.16) related one to another via a composition 
transform (2.17)? 



x = P(t) 



y = Q(t). 



(2.15) 




(2.16) 
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If the answer to this question is positive, we can restrict our parametriza- 
tions P,Q to be "mutually prime in composition sense" (see [12]) and thus 
to obtain uniqueness of the reconstruction. 

More generally, the "inversion problem" for system (2.16) is: 

To characterize all the solutions of system (2.16) and to provide an effec- 
tive way to find these solutions. 

A special case of the inversion problem, in which definite results have been 
recently obtained, is the "Moment vanishing problem" : 

To characterize all the pairs P, Q for which the moments defined by 
(2.16) vanish. 

The moment vanishing problem plays a central role in study of the center 
conditions for the Abel differential equation (see [3j[U[7], [2l]-[36]). In fact, 
it provides an infinitesimal version of the Poincere Center-Focus problem for 
the Abel equation. In spite of a very classical setting (we ask for conditions 
of orthogonality of pQ to all the powers of P!) this problem has been solved 
(for P and Q polynomials) only very recently ([35J). Let us describe the 
solution. 

We say that P and Q satisfy a "composition condition" if there are poly- 
nomials P(w) and Q(w), and a polynomial W(x), satisfying W(0) = W(l), 
such that 

P( x ) = P(W(x)), Q(x) = Q(W(x)). (2.18) 

Composition condition (2.18) can be easily shown to imply the vanishing of all 
the moments (2.16). In many cases it is also a necessary one, but not always. 
The examples of P, Q annihilating the moments (2.16) but not satisfying 
(2.18) can be obtained as follows (see [M]): if P has two right composition 
factors Wi(x) and W2(x), then P and Q = Wi + Wi will annihilate the 
moments (2.16) because of a linearity with respect to Q. For some P we 
can find W\ and W% which are mutually prime in composition algebra (see 
|42j). Then typically P and Q = W% + W2 will have no common right 
composition factors ([M]). The result of |32J claims that this is essentially 
the only possibility: 

Theorem 2.3 (JEM/) All the moments (2.16) vanish if and only if Q is a sum 
°f Qji 3 — 1) • • • ) I: such that for each j the polynomials P and Qj satisfy 
composition condition (2.18). 
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One can expect that the methods developed in [2HHEIE2], can help 

in further analyzing the reconstruction problem for semi-algebraic functions 
in one and more variables. See, in particular, Section 3.3. below. 



3 Functions of two variables 

Also in two dimensions exact reconstruction of semi-algebraic functions (and 
of their extension to SA(tpi, . . . , ipi)) can be considered as one of the ultimate 
goals of our approach. 

3.1 Reconstruction of polygons from complex moments 

In [301 [181 [T3] the problem of reconstruction of a planar polygon from its 
complex moments is considered. The complex moments of a function f(x, y) 
are defined as 



Complex moments can be expressed as certain specific linear combinations 
of the real double moments mki(f)- 

For a plane subset A its complex moments fik(A) are defined by fj,k(A) = 
HkiXA), where \a is the characteristic function of A. 

Let P be a closed n-sided planar polygon with the vertices z%, % — 1, . . . , n. 
The reconstruction method of [30J is based on the following result of [8J: 

Theorem 3.1 There exist a set of n coefficients depending 
only on the vertices Zi, such that for any analytic function <p(z) on P we have 



f>k(f)= / / z k f(x,y)dxdy, k 




0, 1, . . . , z = x + iy. 



(3.1) 




The coefficients a,, j — 1, . . . ,n are given as aj 




z i- z i+i 

Zj-Z j + 1 



Applying this formula to <p(z) = z k we get 



n 



k(k - l)^ k - 2 {xp) = Yl 




0,1,... 



(3.2) 



i=l 
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where we put /i_ 2 = V>-x = 0. So on the left-hand side we have shifted 
moments of P. 

If we ignore the fact that a,- can be expressed through Zi and consider 
both dj and z% as unknowns, we get from (3.2) a system of equations 

n 

Y^a iZ 1 = v k , A; = 0,1,..., (3.3) 
i=i 

where v k denotes the "measurement" k(k — l)/x^_ 2 (P). System (3.3) is iden- 
tical to system (2.7) which appears in reconstruction of linear combination 
of ^-functions. One of the solution methods suggested in [30] is the Prony 
method described in Section 2.4 above. Another approach is based on ma- 
trix pencils. In [18j [13] an important question is investigated of polygon 
reconstruction from noisy data. 



3.2 Quadrature domains 

We introduce, following [19] , a slightly different sequence of double moments: 
for a function g(z) = g(x + iy) the moments rhki(g) are defined by 

m k i(g) = J J z k z l g(z)dxdy, k, I e N. (3.4) 

One defines the moment generating function I g (v,w) = YlTi=o 1 ^ l ki(g)v k w l 
and the "exponential transform" 

I g (v, w) = 1 - exp(—I g (v, w)) = 

TV 




exp( / / - — ' J := V b k i(g)v n w\ 

HI L 17 — 1)11 ■? — in\ 



g(z)dxdy _ x - , , , 
(z — v)(z — w) 

K JK ' k,l=0 

Now (classical) quadrature domains in C are defined as follows: 

Definition 3.1 A quadrature domain Q C C is a bounded domain with 
the property that there exist points Z\, . . . , z m G Q and coefficients c^, i = 
1, . . . ,m, j — 0, . . . , Si — 1, so that for all analytic integrable functions f(z) 
in Q we have 

m s^ — 1 

f(x + iy)dxdy = YsH djfiKzi)- (3-5) 

i=l j=0 

N = s\ + • • • + s m is called the order of the quadrature domain Q. 




16 



The simplest example is provided by the disk D R (0) of radius R centered at 
G C: f f Djt ,Q\ f(x + iy)dxdy = nR 2 f(0). The results of Davis ([8]; Theorem 
3.1 above) give another example in this spirit. 

The following result ([20], [19], Theorem 3.1) provides a necessary and 
sufficient condition for Q C C to be a quadrature domain: let Iq(v,w) = 
I xn {v, w) be the exponential transform of Q. 

Theorem 3.2 Q is a quadrature domain if and only if there exists a poly- 
nomial p(z) with the property that the function q(z,w) = p(z)p(w)In(z,w) 
is a polynomial at infinity (denoted by q(z,w)). In that case, by choosing 
p(z) of minimal degree, the domain Q is given by Q = {z G C, q(z, z) < 0}. 
Moreover, the polynomial p(z) in this case is given by p(z) = YliLi( z ~ z i) Si > 
where Zi are the quadrature nodes ofQ. 

Now, the algorithm in [19] for reconstruction of a quadrature domain from 
its moments consists of the following steps: 

1. Given the moments rh^Q) = rhki(xn) up to a certain order, construct 
the (truncated) exponential transform I(v,w) = X]fc°/=o bkiv k w l . 

2. Identify the minimal integer iV such that det(6/ C) /)^ r /=0 = 0. Then 
there are coefficients ctj, j = 0, . . . , N — 1, such that for B = (b k)[ )^ l=0 and 
a = (ai, . . . , CKjv-i, 1) T we have 

Ba = 0. (3.6) 

We solve this system with respect to a. Then the polynomial p(z) defined 
above is given by p(z) = z N + o.n-iz n ~ 1 + • • ■ + cto- 

3. Construct the function 

1 N ~ l 11 
R n (z,w) = p(z)p(w)exp( ^ w (ft)____) 

k,l=0 

and identify q(z,w) as the part of Rq(z,w) which does not contain negative 
powers of z and w. Then the domain Q is given by Q = {z G C, q(z, z) < 0}. 

Remark. Let us substitute into the definition of the quadrature domain 
(formula (3.5) above) f(z) = z k . Assuming that all the quadrature nodes Zi 
are simple, we get for the complex moments rn.fc j o(^) = rrik{VL) the expression 

rn 

m k {Vt) = ^ CjZj, 
i=i 
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which is identical to (3.3) in reconstruction of planar polygons. So we can 
reconstruct the quadrature nodes Zi and the coefficients q from the complex 
moments only, and we get once more a complex system which is identical to 
(2.7). Allowing quadrature nodes %i of an arbitrary order, we get a system 
corresponding to a linear combination of 5-functions and their derivatives 
(compare [271 129] ). 

Notice that system (3.6) that appears in step 2 of the reconstruction 
algorithm above is very similar to system (2.10) in the solution process of 



3.3 (5-functions along algebraic curves 

As we've mentioned above, a natural class of functions f(x, y) of two vari- 
ables, for which we can hope for an explicit reconstruction from a finite 
number of the moments rriki(f) = J J x k y l f(x,y)dxdy, k, I — 0,1, ... , con- 
sists of semi- algebraic functions. Those are piecewise-algebraic functions 
with the continuity pieces bounded by piecewise-algebraic curves. Among 
semi-algebraic functions are piecewise-polynomials with the continuity pieces 
bounded by spline curves - a very natural and convenient object in construc- 
tive approximation. 

Most of the methods presented in Section 2 for functions of one variable 
are applicable also in the case of two variables. In particular, generalizing 
the approach of [UJ we can differentiate piecewise-polynomials a sufficient 
number of times and finally get a combination of weighted 5 functions along 
the partition curves. See also [29J. 

In this paper we restrict ourself to a discussion of only one example. 
Assume that f(x, y) is a 5-function 5s along a rational curve S, i.e. for any 
ijj(x,y) we have J J fifidxdy = J s tfj(x,y)dx. Let 



be a rational parametrization of S. The moments now can be expressed as 



where p denotes the derivative P' of P. This system is an extension of system 
(2.16): here we are allowed to use all the double moments, while in (2.16) 
only the moments are available. 



(2.7). 



x = P(t), y = Q(t), t e [0,1] 



(3.7) 




(3.8) 
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Also here we notice that a rational parametrization P, Q of the curve S 
in general is not unique: for any rational function W(t) satisfying W(0) = 
0, W(l) = 1 we get another rational parametrization of our curve: 

x = P(t), y = Q{t), with P(t) = P(W(t)), Q(t) = Q(W(t)). (3.9) 

Consequently, we can reiterate the question in Section 2.6 with better chances 
for a positive answer: 

Are all the solutions of (3.9) related one to another via a composition 
transform (3.10)? 

For the "Moment vanishing problem" for (3.9) a definite answer has been 
obtained in [36J: composition condition (2.18) is necessary and sufficient for 
the moments vanishing. 

Let us now assume that the curve S is closed and that it can be parametrized 
by x = P(t), y = Q(t), with t in the unit sircle S 1 . The study of the double 
moments of this form brings us naturally to the recent work of G. Henkin 
[21 1221123]. Indeed, the vanishing condition for the moments (4.12) is given 
by Wermer's theorem ([I]): rriki(f) = if and only if S bounds a complex 
2-chain in C 2 . See [3J for a simple interpretation of Wermer's condition in 
the case of rational P,Q. In general, if the moments rriki(f) do not vanish 
identically, then the local germ of complex analytic curve S generated by S 
in C 2 does not "close up" inside C 2 . G. Henkin's work ( [21 [221 [23] ) , in par- 
ticular, analyzes various possibilities of this sort in terms of the "moments 
generating function" . We expect that a proper interpretation of the results of 
[21 [221 [23] can help also in understanding of the moment inversion problem. 
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